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Acoustic  travel-time  tomography  allows  one  to  reconstruct  temperature  and  wind  velocity  fields  in 
the  atmosphere.  In  a  recently  published  paper  [S.  Vecherin  et  ai,  J.  Acoust.  Soc.  Am.  119,  2579 
(2006)],  a  time-dependent  stochastic  inversion  (TDSI)  was  developed  for  the  reconstruction  of  these 
fields  from  travel  times  of  sound  propagation  between  sources  and  receivers  in  a  tomography  array. 

TDSI  accounts  for  the  correlation  of  temperature  and  wind  velocity  fluctuations  both  in  space  and 
time  and  therefore  yields  more  accurate  reconstruction  of  these  fields  in  comparison  with  algebraic 
techniques  and  regular  stochastic  inversion.  To  use  TDSI,  one  needs  to  estimate  spatial-temporal 
covariance  functions  of  temperature  and  wind  velocity  fluctuations.  In  this  paper,  these 
spatial-temporal  covariance  functions  are  derived  for  locally  frozen  turbulence  which  is  a  more 
general  concept  than  a  widely  used  hypothesis  of  frozen  turbulence.  The  developed  theory  is  applied 
to  reconstruction  of  temperature  and  wind  velocity  fields  in  the  acoustic  tomography  experiment 
carried  out  by  University  of  Leipzig,  Germany.  The  reconstructed  temperature  and  velocity  fields  are 
presented  and  errors  in  reconstruction  of  these  fields  are  studied.  ©  2007  Acoustical  Society  of 
America.  [DOI:  10.1121/1.2756798] 
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I.  INTRODUCTION 

Knowledge  about  temperature  and  wind  velocity  fields 
is  important  in  many  disciplines,  e.g.  boundary  layer  meteo¬ 
rology,  theories  of  turbulence,  and  studies  of  electromagnetic 
and  acoustic  wave  propagation  in  the  atmosphere.  A  conven¬ 
tional  way  to  measure  the  temperature  and  wind  velocity  is 
to  use  in  situ  sensors.  However,  such  measurements  can  pro¬ 
vide  the  data  only  at  some  spatial  points.  To  measure  the 
temperature  and  velocity  fields  with  high  resolution,  a  very 
large  number  of  sensors  is  required.  This  will  not  only  in¬ 
crease  the  cost  of  measurements  but  also  distort  the  original 
fields.  This  problem  can  be  overcome  by  using  acoustic  to¬ 
mography. 

Acoustic  travel-time  tomography  of  the  atmosphere  al¬ 
lows  one  to  reconstruct  temperature  T( R)  and  wind  velocity 
V(R)  fields  in  a  tomographic  volume  or  area  given  positions 
of  sound  sources  and  receivers  and  travel  times  f'r  of  sound 
propagation  between  them.  Here,  R  =  (x,y,z)  is  a  vector  in 


the  Cartesian  coordinate  system  and  1=1,2,...,/,  where  I  is 
the  number  of  data  tf,  i.e.,  the  number  of  sound  propagation 
paths  between  sources  and  receivers. 

In  the  past  decade,  several  acoustic  travel-time  tomog¬ 
raphy  experiments  have  been  carried  out  in  the  lowermost 
few  meters  of  the  atmosphere.1  9  Further,  some  works  were 
done  on  numerical  simulation  of  acoustic  tomography  of  the 
atmosphere. llW  l  One  of  the  main  problems  in  both  tomogra¬ 
phy  experiments  and  numerical  simulations  is  to  find  a  good 
inverse  algorithm  for  reconstruction  of  T  and  V  fields.  Most 
of  the  algorithms  used  so  far  employed  a  partition  of  a  to¬ 
mographic  volume  into  grid  cells  where  the  values  of  T  and 
V  are  constant,  and  a  subsequent  solution  of  a  set  of  alge¬ 
braic  equations.3^9'11’12’14’15  Numerical  simulations  have 
shown  “  that  these  algebraic  algorithms  can  give  a  good  re¬ 
construction  of  7  and  Y  fields  if  the  number  of  the  grid  cells 
is  small  enough  so  that  a  corresponding  inverse  problem  is 
overdetermined  (the  number  of  equations  is  greater  than  the 
number  of  unknowns).  However,  in  two-dimensional  (2D) 
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tomography  experiments,  for  example,  there  are  three  con¬ 
tinuous  fields  to  reconstruct  (temperature  and  two  compo¬ 
nents  of  velocity),  whereas  the  number  of  data  is  limited.  (A 
large  number  of  sources  and  receivers  can  distort  the  T  and  Y 
fields  inside  a  tomographic  area.)  Therefore,  the  number  of 
grid  cells  where  T  and  V  are  reconstructed  is  relatively  small 
(a  few  dozen  or  less). 

To  increase  spatial  resolution  of  the  reconstructed  tem¬ 
perature  T  and  velocity  V  fields,  one  can  use  a  stochastic 
inversion  (SI),  which  is  appropriate  for  solutions  of  inverse 
problems  involving  continuous  fields.2'15  In  SI,  the  increase 
in  spatial  resolution  comes  with  a  price:  One  needs  to  know 
spatial  covariance  functions  B^R^  ,R2)  and  fiI-(R1,R2)  of 
temperature  and  velocity  fluctuations  at  the  time  moment  t  of 
an  experiment.  Here,  s  stands  for  spatial,  and  i  and  j  indicate 
the  components  of  wind  velocity  vector.  Recent  studies  have 
shown  ’  that,  in  SI,  the  use  of  different  covariance  func¬ 
tions  (exponential,  Gaussian,  von  Karman,  and  some  others) 
results  in  the  reconstructed  temperature  and  velocity  fields 
which  are  nearly  identical  (with  a  relationship  between  cor¬ 
relation  lengths  that  provides  the  equality  of  the  integral 
lengths  of  these  functions,  see  Ref.  18  for  details).  Therefore, 
one  of  these  covariance  functions,  e.g.,  Gaussian,  can  be 
used  in  SI.  Note  that  in  SI  the  accuracy  of  T  and  Y  recon¬ 
struction  increases  with  the  increase  of  the  number  I  of  data. 

To  increase  the  number  of  data  without  increasing  the 
number  of  sources  and  receivers,  a  time-dependent  stochastic 
inversion  (TDSI)  in  acoustic  tomography  of  the  atmosphere 
was  proposed  in  Ref.  19.  In  TDSI,  the  travel  times  ff  are 
measured  repeatedly  at  the  time  moments  t\,t 2, ... , tn.  Then, 
the  fields  T(R,t)  and  V( R ,  t)  are  reconstructed  given  n  sets 
of  tf,  positions  of  sources  and  receivers,  and  assuming  that 
the  spatial-temporal  covariance  functions  of  temperature  and 
wind  velocity  fluctuations,  /i77(R|  ,  q ;  R2,i2)  and 
^//(R  i  Ti  ;R2T2),  are  known.  In  Ref.  19,  these  spatial- 
temporal  covariance  functions  were  determined  assuming 
that  turbulence  is  frozen,"  i.e.,  temperature  and  wind  veloc¬ 
ity  fluctuations  are  advected  with  the  constant  wind  velocity 
V0.  Numerical  simulation  of  acoustic  tomography  of  the  at¬ 
mosphere  showed  that  TDSI  allows  better  reconstruction  of 
temperature  and  velocity  fields  than  SI  does.19 

A  general  idea  to  use  spatial-temporal  covariances  in  SI 
is  known  in  the  literature.  For  example,  this  idea  has  been 
successfully  used  in  the  satellite  altimetry  of  the  ocean  sur- 
face  level."  TDSI  in  acoustic  tomography  of  the  atmo¬ 
sphere  is  somewhat  similar  to  methods  used  in  the  satellite 
altimetry.  However,  a  mathematical  apparatus  of  TDSI  is  dif¬ 
ferent  from  that  used  in  altimetry  and  was  developed  in  Ref. 
19. 

The  main  goal  of  the  present  paper  is  further  develop¬ 
ment  of  a  mathematical  apparatus  of  TDSI  and  the  use  of  the 
developed  theory  in  the  reconstruction  of  the  T  and  V  fields 
in  the  acoustic  tomography  experiment  STINHO  (STructure 
of  turbulent  transport  under  INHOmogeneous  surface  condi¬ 
tions)  carried  out  by  University  of  Leipzig,  Germany.9  First, 
we  derive  analytical  formulas  for  B^Rj  ,t1  ;R2,r2)  and 
^i/(R  i  Ti  ;R2T2)  for  the  case  of  locally  frozen  turbulence 
which  is  a  generalization  of  the  hypothesis  of  frozen  turbu¬ 


lence.  The  obtained  formulas  for  BTT(Rl,t1',R2,t2)  and 
Bjj(Rl  Ji  ;R2,f2)  account  for  the  variance  of  wind  velocity 
fluctuations  crv  and,  hence,  are  much  more  realistic  than 
those  for  the  case  of  frozen  turbulence  where  oy  is  tacitly 
assumed  0.  Second,  we  determine  errors  in  the  input  data  for 
TDSI.  Finally  using  these  errors  and  formulas  for 
Btt^R]  , z i  ;R2,f2)  and  Bjj  (R!,ri;R2,r2)  for  locally  frozen  tur¬ 
bulence,  we  apply  TDSI  for  reconstruction  of  T  and  V  fields 
in  the  acoustic  tomography  experiment  STINHO.  Note  that 
this  is  the  first  reconstruction  of  the  temperature  and  wind 
velocity  fields  in  an  acoustic  tomography  of  the  atmosphere 
experiment  using  TDSI. 

The  paper  is  organized  as  follows.  Section  II  describes 
starting  equations  for  acoustic  tomography  and  the  TDSI  al¬ 
gorithm.  Formulas  for  spatial-temporal  covariance  functions 
of  temperature  and  wind  velocity  fluctuations  for  locally  fro¬ 
zen  turbulence  are  derived  in  Sec.  III.  Using  the  residual 
analysis,  the  errors  in  the  input  data  for  TDSI  are  obtained  in 
Sec.  IV.  In  Sec.  V,  new  methodology  is  applied  to  reconstruc¬ 
tion  of  temperature  and  velocity  fields  in  the  acoustic  tomog¬ 
raphy  experiment  STINHO.  The  conclusions  are  presented  in 
Sec.  VI. 


II.  STARTING  EQUATIONS 

In  this  section,  starting  equations  for  acoustic  travel-time 
tomography  of  the  atmosphere  and  a  brief  description  of 
TDSI  are  presented.  A  detailed  description  of  TDSI  can  be 
found  in  Ref.  19. 

Let  u{R,t),  u(R,/),  and  vv( R , t)  be  x,  y,  and  z  compo¬ 
nents  of  the  random,  three-dimensional  (3D)  vector  of  wind 
velocity  V(R,f): 

V(R,f)  =  i7(R,r)et  +  u(R,r)ev  +  w(R,t)e,.  (1) 

Here  et,  ev,  and  e,  are  the  unit  vectors  along  the  x,  y,  and  z 
axes,  respectively.  The  adiabatic  sound  speed  cL,  temperature 
T,  and  wind  velocity  V  within  a  tomographic  area  at  the  time 
moment  t  can  be  represented  as  sums  of  their  spatially  aver¬ 
aged  values  c0(t),  T0(t),  Y0(t)  =  {u0(t)  ,v0(t)  ,w0(t))  and 
their  fluctuations  c(R,t),7’(R,t),V(R,f)  =  (M(R,/),u(R,t), 
w(R,f)): 

G.(Rt)  =  C0(f)  +  c(R,t),T  =  T0(t)  +  T(R,t), 

(2) 

V(R,f)  =  V0(r)  +  V(R,r). 

Note  that  tomography  allows  one  to  reconstruct  acoustic  vir¬ 
tual  temperature,  which  we  denote  as  T.  It  is  related  to  the 
thermodynamic  temperature  Tt  h  by  the  following  relation¬ 
ship:  T=Tth(  1+0  .51  Iq),  where  q  is  the  specific  humidity  of 
■  18 
air. 

We  will  assume  that  the  spatially  averaged  values  are 
equal  to  the  mathematical  expectations  of  the  random  fields 
at  each  spatial  point  at  the  time  moment  f,  i.e., 
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-  COS  <Pi 


-  sin  c pj 


(8) 


(cL(R,t))  =  c0(t),  (T(R,t))  =  T0(t),  <V(R,f)>  =  V0(f). 

(3) 

Hereinafter,  the  angular  brackets  ( )  denote  averaging  over  an 
ensemble  of  realizations.  Note  that  it  follows  from  Eqs.  (2) 
and  (3)  that  the  mathematical  expectations  of  the  fluctuations 
at  any  point  R  and  time  t  are  equal  to  zero. 

Most  acoustic  tomography  arrays  used  so  far  were  2D. 
In  Ref.  19  TDSI  was  developed  for  the  2D  case  in  which 
sources,  receivers,  and  sound  propagation  paths  are  located 
in  one  plane.  For  concreteness,  we  will  assume  that  this 
plane  coincides  with  the  (x ,  y)  plane.  In  this  case,  the  travel 
times  of  sound  propagation  between  sources  and  receivers 
are  given  by  linearized  equations  of  the  forward  problem  of 
acoustic  travel-time  tomography  of  the  atmosphere19: 


<T«  = 


h 

i 


1  - 


Uo(t)sxi  +  v0  (t)svl 


c0(t ) 


clit) 


HI 


c0(t) 
2  T0(t) 


T{R,t)  +  u(R,t)sxl  +  v{R,t)syi  f . 


(4) 


Here  L,  is  the  distance  between  a  source  and  receiver,  sxl 
=  cos  (ph  and  sv/=sin  <ph  where  tpt  is  the  angle  between  ex  and 
the  direction  of  sound  propagation.  Note  that  in  the  linear¬ 
ized  forward  problem  the  paths  of  sound  impulses  can  be 
approximated  by  straight  lines. 

The  inverse  problem  of  2D  acoustic  travel-time  tomog¬ 
raphy  of  the  atmosphere  is  to  reconstruct  the  mean  fields 
T0(t),  u0(t),  Void  and  their  fluctuations  T(R,t),  u(R,t), 
u(R,f)  given  the  travel  times  tf{t),  the  angles  <pit  and  the 
distances  Lr  This  reconstruction  consists  of  two  steps.  First, 
Tg{t),  ug{t),  and  Vp(t)  are  reconstructed  using  the  least- 
squares  estimation.  Then,  the  fluctuations  T(R,t),  u{R,t), 
and  u(R,r)  are  reconstructed  with  the  help  of  TDSI. 

When  reconstructing  T0(t),  n0(r),  and  Up(t),  the  integral 
term  on  the  right-hand  side  of  Eq.  (4)  is  neglected.  In  this 
case,  Eq.  (4)  can  be  rewritten  in  a  matrix  notation: 


Gf  =  b. 


(5) 


Here  the  elements  of  the  column-vector  b  are  known  (can  be 
calculated  from  experimentally  measured  values  of  tf(t)  and 

Li)- 


(6) 


the  unknown  column-vector  f  has  three  elements 


1  _  up (t) 

Cp(t )  ’  '  2  Co(f)  ’ 


and  the  matrix  G  is  given  by 


(7) 


G  = 


1  -  COS  (Pj 


-  sin  tpi 


Then,  the  overdetermined  inverse  problem  for  the  elements 
of  vector  f  (/  is  assumed  to  be  greater  than  3)  is  solved  using 
the  least-squares  estimation: 


f=(GrG)_1Grb. 


(9) 


Here  f  is  the  optimal  least-squares  estimation  of  f.  The  opti¬ 
mal  estimations  c0,  u0,  v0  of  the  mean  fields  c0,  up,  Vp  can 
easily  be  calculated  from  f  using  Eq.  (7).  The  mean  tempera¬ 
ture  T0  can  be  obtained  from  c0  using  the  following  formula: 

cl=yRaTp,  (10) 


where  y~  1.41  is  the  ratio  of  the  specific  heats  and  Ra  is  the 
gas  constant  for  dry  air. 

The  estimated  values  Cp,  T0,  u0,  and  v0  are  substituted 
back  into  Eq.  (4).  The  resulting  equation  can  be  written  in 
the  following  form: 


<7,0)  =  J 


dl{ 


£p(t ) 

2  T0{t) 


r(R,f)  +  w(R,t)ij:(  +  i;(R,f)s  ■  ?■ , 


where  q^i)  are  given  by 

qM  =  L[c0(t )  -  u0(t)sxi  -  Vp(t)iyi]  -  cl{t)t’[{t) . 


(ID 

(12) 


After  the  reconstruction  of  the  mean  fields,  <7, -(f)  are  known. 

Equation  (11)  is  a  starting  equation  for  TDSI  whose  goal 
is  to  reconstruct  T(R,t),  m(R , r)  and  HR, t)  given  the  values 
of  £/,(r),  <Pr  an(l  Lf.  The  <7, -(f)  are  known  at  the  time  moments 
ti ,  t2 ,  ■  ■  ■ ,  tn  of  an  experiment  and  are  used  to  form  the  vector 
d  of  the  input  data  for  TDSI: 


d  =  [q(f1);q(f2);  ...  ;q(f„)]. 


(13) 


Here  q(4)  are  column  vectors  with  elements  given  by  Eq. 
(12),  and  the  semicolons  indicate  that  these  q  are  arranged  in 
the  column  vector  d. 

Let  m(f)  be  a  column  vector  of  the  fluctuations  being 
reconstructed  at  certain  spatial  points  R^  ...  ,Ry,  where  J  is 
a  total  number  of  points: 

m(f)  =  [T(Rut)  -  ...  ;r(Ry,f);M(R1,r);  ...  ;k(R/T); 

v(Rht)-,...-,v(Rj,t)l  (14) 


Note  that  in  TDSI  the  fluctuations  can  be  reconstructed  at 
any  time  moment  t.  Of  course,  one  can  expect  a  good  recon¬ 
struction  of  T(R,f),  u(R,t),  and  u(R,f)  fields  if  t  is  within 
the  time  interval  [fi,t„]  or  close  to  it.  In  TDSI,  the  optimal 
reconstruction  m(r)  of  m(r)  is  given  by  the  following  for¬ 
mula: 

™(f)  =  CmdCjjd,  (15) 

where  Cmd=(md7)  is  a  model-data  covariance  matrix  and 
C(ld  =  (dd7)  is  a  data  covariance  matrix  given  by  Eqs. 
(18)— (21)  from  Ref.  19.  The  elements  of  these  matrices  are 
expressed  in  terms  of  the  integrals  along  the  sound  propaga- 
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tion  paths  between  sources  and  receivers,  whose  integrands 
contain  the  spatial-temporal  covariance  functions  of  tem¬ 
perature  and  wind  velocity  fluctuations,  Bj-^R]  ,tl  ;R2,f2) 
and  Bjj( Rj  ,tl ;  R2 , t2) .  In  an  experiment,  the  components  q,  of 
the  vector  d  are  always  known  with  some  uncertainties  (ob¬ 
servation  noise)  which  we  denote  as  Aoi.  As  a  result,  the 

'  IQ 

matrix  Cdd  has  additional  terms  on  its  main  diagonal  equal 
to  the  variances  of  these  uncertainties  ct^=(A4(.). 

Expected  squared  errors  of  the  estimation  m(r)  are  ele¬ 
ments  on  the  main  diagonal  of  the  error  covariance  matrix 
r,  .  ,19 

Cee,  given  by 

c«=cmm-cmdcd‘cL,  (16) 


T(R,t2)  =  T(R-Y(R,t1)T,h),  (20) 

where  T=t2-tx.  Note  that  in  Eq.  (20)  the  velocity  Y  is  a 
random  function  of  R  and  f, .  If  in  Eq.  (20)  Y  were  constant, 
this  equation  would  become  a  formula  for  frozen  turbulence 
which  has  been  widely  used  in  the  literature."  Thus,  the 
hypothesis  of  locally  frozen  is  a  generalization  of  that  of 
frozen  turbulence. 

Equation  (20)  allows  one  to  express  the  spatial-temporal 
covariance  function  Bj-^R] , f x :  R2,f2)  in  terms  of  the  spatial 
covariance  function  B^T(R1,R2).  Using  this  equation,  we 
have: 


where  Cmm=(mm7)  is  a  model  covariance  matrix. 

Thus,  in  order  to  implement  TDSI  for  reconstruction  of 
temperature  and  wind  velocity  fields,  one  needs  to  know 
spatial-temporal  covariance  functions  Bj^R^y  ;R2,f2)  and 
B^Rj  ,tx  ;R2,r2)  and  the  variances  cr(.  These  covariance 
functions  and  variance  will  be  calculated  in  the  next  two 
sections. 

III.  LOCALLY  FROZEN  TURBULENCE 

In  this  section,  analytical  formulas  for  the  spatial- 
temporal  covariance  functions  of  temperature  and  wind  ve¬ 
locity  fluctuations,  B^R]  ,t\  ;R2,t2)  and  B,;(R|  ;R2,f2), 

are  derived  for  the  case  of  locally  frozen  turbulence. 

For  this  derivation,  we  will  need  particular  forms  of  the 
spatial  covariance  functions  of  temperature  and  wind  veloc¬ 
ity  fluctuations  B'^Rj  ,R2)  and  B^(R1;R2).  We  will  assume 
that  turbulence  is  statistically  isotropic  and  homogeneous, 
and  that  the  spatial  covariance  functions  of  temperature  and 
longitudinal  velocity  fluctuations  are  given  by  Gaussian 
functions.  The  components  of  the  tensor  BF(R1;R2)  can  be 
calculated  from  the  covariance  function  of  longitudinal  ve¬ 
locity  fluctuations  using  Eqs.  (6.34)  and  (6.35)  in  Ref.  18.  As 
a  result,  we  obtain  the  following  spatial  covariance  functions 
of  temperature  and  velocity  fluctuations: 

Btt(Ri,R2)  =  crj-exp^-  j ,  (17) 

B^Rt.Rz)  =  4  exP(-  j}  ( 1  -  ,  (18) 

B^(R1,R2)  =  4exp(-|j^.  (19) 

Here,  oy  and  crv  are  the  standard  deviations  of  temperature 
and  velocity  fluctuations,  lT  and  lv  are  their  correlation 
lengths,  and  p=R2-R1  =  (pt,pv,p.).  Other  components  of 
the  tensor  B^  (R,,R2)  are  given  by  formulas  similar  to  those 
for  B'h(R,,R2)  and  B*u(R],R2). 

A  hypothesis  of  locally  frozen  turbulence  is  formulated 
in  Ref.  24.  According  to  this  hypothesis,  during  a  relatively 
small  time  interval  r,  velocities  V(R .  t)  can  be  considered  as 
constant  in  the  vicinity  of  each  spatial  point.  This  allows  one 
to  express  the  temperature  field  T( R ,  t2)  at  the  time  moment 
t2  in  terms  of  the  temperature  field  at  the  time  moment  t\ : 


Br/(R  1  ,f  1 ;  R2,f2)  —  (T(Ri,ti)7f  RiT^^y 

=  (T(Ri,ti)T(R2  -  V(R2,f1)r,f1))jjv. 

(21) 

Here,  ( )Ty  denotes  the  average  over  an  ensemble  of  realiza¬ 
tions  of  the  random  fields  T  and  V.  These  fields  are  uncor¬ 
related  in  the  model  of  homogeneous  isotropic  turbulence 
which  is  used  in  this  paper.-4  Assuming  Gaussian  probability 
density  distribution,  these  random  fields  are  statistically  in¬ 
dependent.  Therefore,  in  the  right-hand  side  of  Eq.  (21),  one 
can  first  calculate  the  average  over  an  ensemble  of  realiza¬ 
tions  of  T  which  yields  the  spatial  covariance  function  B'rl 
by  definition.  Taking  into  account  that  p=R2-R],  we  have 

Br7(p,T)  =  <BJrr(p-V(R2TiK)>v-  (22) 


The  average  on  the  right-hand  side  of  this  formula  is  calcu¬ 
lated  in  the  Appendix.  The  result  is 


Br7-(p,  t)  =  <7^.  exp 


(P~  Vq(tK)2 

?T 


(23) 


where  the  effective  variance  crj.  and  the  square  of  correlation 


length  lj  are  given  by 


cr^  = 


1+2 


3/2’ 


Tj=12t+2o1vtl. 


(24) 


Note  that  Eqs.  (23)  and  (24)  can  also  be  derived  from  Eq. 
(31.29)  in  Ref.  24.  The  latter  equation  gives  a  formula  for  a 
Fourier  transform  of  Brr(p,  r)  with  respect  to  p  and  r,  valid 
for  an  arbitrary  spatial  covariance  function  and  probability 
density  distribution. 

Similarly  to  the  derivation  of  Eq.  (23),  one  can  obtain  a 
formula  for  the  spatial-temporal  covariance  function  of  the 
longitudinal  velocity  fluctuations  and,  then,  formulas  for 
(p,r)  and  Buv(pj). 

(P~  Vq(T)t')2 

lv  J 

x  /  !  _  (py  ~  t/||(t|  )  ~)~  +  (p:  -  U’o(t|  )  T)2 

\  Fy 

(25) 


B,Jp,t)  =  a2vexp 
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(p-  v0(?i)t)2 


(31) 


=  of,  exp 


lv 


X 


(px-  Uoit^T^py-Voitijr) 


f- 

lv 


where 


ai 


1+2 


3/2’ 


?v=l2v+2a2vT2. 


(26) 


(27) 


Equations  (23)-(27)  provide  analytical  formulas  for  the 
spatial-temporal  covariance  functions  of  temperature  and  ve¬ 
locity  fluctuations  for  the  considered  case  of  locally  frozen 
turbulence.  In  the  limiting  case  crv=  0  and  time-independent 
V0,  these  formulas  coincide  with  those  for  frozen 
turbulence.19  It  follows  from  Eqs.  (23)— (27)  that  the  depen¬ 
dence  of  the  spatial-temporal  covariance  functions  on  r 
manifests  in  three  effects:  The  spatial  arguments  of  the  co- 
variance  functions  are  shifted  by  the  vector  V0(f1)r,  and  the 
effective  variances  of  the  fluctuations  decrease  while  the  ef¬ 
fective  correlation  lengths  increase  with  the  increase  in  r. 

A  necessary  condition  of  applicability  of  the  hypothesis 
of  locally  frozen  turbulence  is  formulated  in  Ref.  24: 

aV«|V0|.  (28) 


In  this  paper,  we  consider  the  case  when  the  mean  wind 
velocity  V0  can  slowly  depend  on  time  t.  This  imposes  an 
additional  necessary  condition  of  applicability  of  the  hypoth¬ 
esis  of  locally  frozen  turbulence: 

|VoW-V0(f1)|<|V0(fi)|,  (29) 

where  t  e  [/, ,  f2]-  As  far  as  we  know  sufficient  conditions  of 
applicability  of  the  hypothesis  of  locally  frozen  turbulence 
have  not  been  obtained.  They  should  probably  include  a  limi¬ 
tation  on  the  time  interval  r=t2-t t  during  which  this  hypoth¬ 
esis  is  valid. 


IV.  ESTIMATION  OF  THE  ERRORS 

In  this  section,  the  variances  a 2qi  of  the  uncertainties  Af/ 
in  the  input  data  d  are  calculated.  It  follows  from  Eq.  (11) 
that  these  variances  depend  on  the  variances  of  mean  squared 
errors  cif,  a2/,  of,  and  a2  in  reconstruction  of  c0,  70,  u0,  and 
v0.  Therefore,  we  first  estimate  these  errors. 

Let  si  be  the  estimated  variance  of  components  b,  of  the 
vector  b  which  appears  in  the  right-hand  side  of  Eq.  (5). 
Using  a  residual  estimation,16  we  have 

2  lb  -  Gfl  'fb  -  Gf  )  ,  . 

4  = - JZj - Z>  00) 

where  1-3  is  a  number  of  independent  degrees  of  freedom  in 

2 

the  numerator  of  Eq.  (30).  The  mean  squared  errors  cr~  of  the 
vector  f  are  diagonal  elements  of  its  covariance  matrix16  cr~ 
=  diag[Rjr]-  Here 


r;.;.=^(g/G)-1. 


Using  Eqs.  (7)  and  (10)  and  neglecting  the  terms  of  order 
(uq/cq)2  and  ( Vq/cq )2,  the  mean  squared  errors  of ,  erf,  cif, 
and  erf  can  be  expressed  in  terms  of  the  mean  squared  errors 


d2-- 


C°°/V 


6\- 


4  c2 


iyRa) 


*4  2 
c°°/2’ 


of: 


.4  2 

C0°/3- 


(32) 


Once  the  errors  cif,  erf,  erf,  and  of  are  known,  it  is  pos¬ 
sible  to  estimate  the  variances  <r2  of  the  errors  in  the  input 
data  qj(t)  for  TDSI.  The  errors  in  q,(t)  are  due  to  errors  in  Lh 
sxi,  and  syi  caused  by  uncertainties  in  transducers  positions, 
errors  in  measurements  of  t1',  and  errors  in  the  reconstruction 
of  c0(t),  f0(t),  u0(t ),  and  v0(t).  As  L?  =  (xi-x0i)2  +  (yi-y0i) 2, 
where  (xhy,)  and  (x0i,y0l)  are  the  coordinates  of  the  receiver 
and  source,  respectively,  one  can  calculate  the  uncertainties 
in  the  sound  propagation  paths  AL  taken  into  account  that 
cos  cpi=(x,-x0i)/Ll  and  sin  <p,-=(y,— yod/Lf. 


;[(A, 


■  Ax0,.)cos  <p,  +  (Av 


'  A'’o,)sin  ft]- 


(33) 


Here,  Av,  Ar  .  A,,,  and  A,,  are  uncertainties  in  the  trans- 

’  xi  x0i  -vi  -v0  i 

ducer  positions.  Therefore,  from  Eq.  (12)  the  uncertainty  Aq 
is 


\  =  (A.tj  -  Aa0,.)(co  cos  (fj  -  u0)  +  (Ay.  -  Ayo.) 

x  (co  Sin  ip,  -  v0)  +  ( Li  -  2c0?f)AC() 

-  Li(\  cos  (pi  +  A„o  sin  (pi)  -  CqA,.,  (34) 

and,  as  m0<Cc0,  v0  <Cc0,  an£l  c0tf~Lh  the  variance 
=  (A^}  of  the  errors  in  q ,•  is 

Ogt  “  2ct^c5  +  Lr(d^  +  dru  cos2  < Pi  +  a 2V  sin2  <p,)  +  c^a2, 

(35) 

where  cTr  is  the  variance  of  the  errors  in  the  transducer  posi¬ 
tions  (ct2  =  (A2)  =  (A2  )  =  (A2 )  =  (A2  )),  and  of  is  the  variance 

r  xj  xqj  yJ  yoi  1 

of  an  error  in  measurement  of  r- . 

V.  ACOUSTIC  TOMOGRAPHY  EXPERIMENT 

In  this  section,  a  developed  theory  of  TDSI  is  applied  to 
reconstruction  of  temperature  T  and  wind  velocity  V  fields  in 
the  acoustic  tomography  experiment  STINHO  carried  out  by 
scientists  from  the  University  of  Leipzig.9  The  experiment 
was  a  part  of  a  larger  meteorological  experiment  to  study 
turbulence,  turbulent  fluxes,  and  other  meteorological  param¬ 
eters  over  heterogeneous  surface. 

The  tomography  array  consisted  of  eight  sources 
SI , . . .  ,S8  and  twelve  receives  R1 , ...  ,R12,  see  Fig.  1.  The 
size  of  the  array  was  300  m  X  440  m.  The  sources  and  re¬ 
ceivers  were  located  above  the  ground  at  an  average  height 
of  2  m.  Most  of  the  ground  within  the  tomography  array  was 
covered  with  grass  except  for  a  spot  of  bare  soil  which  was 
in  the  lower  left  corner  of  the  array  with  the  size  of  90  m 
X  300  m.  This  inhomogeneity  might  create  the  fluctuation  of 
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FIG.  1.  (Color  online)  The  layout  of  sources  (S1....,S8)  and  receivers 
(R1 _ _ R12)  in  the  acoustic  tomography  experiment  STINHO  (Ref.  9). 

temperature  and  wind  velocity  fields  which  is  subject  of  re¬ 
construction  by  TDSI.  Besides,  the  difference  between  mean 
temperatures  over  different  surfaces  was  negligible  (unlike 
the  difference  in  turbulent  heat  fluxes),  and  horizontal  wind 
velocities  were  quite  close  (see  Fig.  2  in  the  Ref.  9).  The 
travel  times  t'f  of  sound  propagation  between  sources  and 
receivers  were  measured  every  minute  during  6  July  2002. 


Our  reconstruction  of  temperature  T  and  wind  velocity 
V  fields  in  this  acoustic  tomography  experiment  consisted  in 
reconstruction  of  their  mean  values  T0(t),  uo(t),  and  u0(t)  and 
fluctuations  T(R,t),  u(R,t),  and  u(R,f)  as  described  in  Sec. 
II.  The  input  experimental  data  for  reconstruction  of  f0(t), 
u0(t),  and  L>n(r)  were  tj,  Lh  and  <p;.  The  reconstructed  mean 
values  of  temperature  and  two  components  of  wind  velocity 
for  the  10  min  time  interval  (5:26-5:35  a.m..  Coordinated 
Universal  Time  (UTC),  which  is  also  known  as  Greenwich 
Mean  Time)  are  plotted  in  Fig.  2  as  solid  lines.  Using  Eq. 
(32),  the  errors  ay,  irlt.  and  &v  in  reconstruction  of  the  mean 
fields  were  calculated.  These  errors  are  plotted  in  Fig.  2  as 
vertical  bars. 

Then,  the  obtained  values  of  T0(t),  «0(r),  and  u0(f)  and 
experimental  values  of  f)1  were  used  to  calculate  the  compo¬ 
nents  <7,(f)  of  the  input  data  d  for  reconstruction  of  fluctua¬ 
tions  7XR,f),  u(R,t),  and  v ( R , t)  using  TDSI.  To  apply 
TDSI,  we  also  needed  to  know  the  variances  crq  of  errors  in 
<7, -(f),  and  parameters  oy,  oy,  lT,  and  lv  of  the  covariance 
functions  BTT(p ,  r)  and  B,j(p,  t).  The  variances  <r2  were  cal¬ 
culated  using  Eq.  (35),  where  dy  dy  and  dr  were  obtained 
from  Eq.  (32),  whereas  a2  and  crj  were  estimated  in  the 
tomography  experiment9  to  be  as  follows:  cr,.=3  cm  and  cr, 
=  0.3  ms.  The  variances  of  temperature  and  velocity  fluctua¬ 
tions,  oy  and  oy,  were  assumed  to  be  the  same  as  the  vari¬ 
ances  measured  in  situ  by  a  sonic  thermometer-anemometer: 
oy=0.27°C  and  oy=0.28  m/s.  The  correlation  lengths  lT 
and  lv  were  chosen  as  follows:  lT=lv=  75  m.  The  values  of  IT 
and  lv  determine  a  characteristic  size  of  fluctuations  7  (  R .  f ) , 
d(R,f),  and  d(R,f)  which  can  be  resolved  in  reconstruction. 
Note  that,  for  the  considered  experiment,  the  values  of  lT  and 
ly  could  not  be  chosen  significantly  less  than  75  m  since  the 
errors  in  reconstruction  of  the  fluctuations  increase  with  de¬ 
creasing  lT  and  lv.u  Finally,  to  apply  TDSI  to  the  reconstruc¬ 
tion  of  the  fluctuations  T(R,t),  d(R,f),  d(R,f),  three  con- 


15.8  - 
15.6  L 


5:26  5:27  5:28  5:29  5:30  5:31  5:32  5:33  5:34  5:35 


FIG.  2.  (Color  online)  The  reconstruc¬ 
tion  of  the  mean  fields  at  5:26-5:30 
a.m.  UTC  on  6  July  2002  by  TDSI 
(solid  lines)  and  SIRT  (dashed  lines). 
The  vertical  bars  indicate  the  esti¬ 
mated  errors  of  mean  field  reconstruc¬ 
tion  by  TDSI. 
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FIG.  3.  (a)  The  TDSI  reconstruction  of  temperature  (°C)  and  wind  velocity  (m/s)  fields  at  5:30  a.m.  UTC  on  6  July  2002.  (b)  The  10  min  averaged  fields  of 
temperature  fluctuations  (°C)  reconstructed  by  (left)  TDSI  and  (right)  SIRT.  The  black  lines  in  the  TDSI  reconstruction  show  the  direction  of  reconstructed 
mean  wind  during  the  interval  of  averaging. 


secutive  sets  of  the  travel  times  tf  were  used.  (In  future 
tomography  experiments,  it  would  be  desirable  to  increase 
the  number  of  such  sets  by  decreasing  a  time  interval  be¬ 
tween  measurements  of  tf) .  After  the  fluctuations  T(  R ,  f) , 
u(R,t),  and  u  (R , t)  were  reconstructed,  the  total  fields  T,  u, 
and  v  were  obtained  as  the  sums  of  the  mean  fields  and 
fluctuations. 

Figure  3(a)  shows  the  fields  T,  u,  and  v  reconstructed  at 
5:30  a.m.  UTC  (In  the  reconstruction,  the  travel  times  tf 
measured  at  5:29,  5:30,  and  5:31  a.m.  were  used.)  “Cold” 
and  “warm”  temperature  eddies  with  different  scales  can  be 
seen  in  Fig.  3(a).  High  and  low  speed  velocity  eddies  can 
also  be  seen.  Figure  4  contains  the  expected  root  mean 
squared  errors  (RMSE)  of  the  reconstruction  of  temperature 
and  wind  velocity  fields  shown  in  Fig.  3(a).  These  errors  are 
calculated  with  the  use  of  Eq.  (16).  The  spatially  averaged 
RMSE  of  temperature  reconstruction  is  0.36  °C.  This  error  is 
less  than  the  temperature  differences  between  cold  and  warm 
eddies  which  can  be  seen  in  Fig.  3(a).  Therefore,  these  cold 
and  warm  eddies  are  reliably  reconstructed.  The  spatially 
averaged  RMSE  of  u  and  v  reconstruction  are  0.35  and 
0.25  m/s,  respectively.  These  errors  are  also  less  than  the 
velocity  differences  between  the  low  and  high  speed  eddies 
which  can  be  seen  in  Fig.  3(a). 

The  reconstructed  values  of  temperature  were  compared 
with  those  measured  by  two  in-situ  sensors  located  within 
the  tomographic  area.  The  results  of  this  comparison  and  the 
coordinates  of  the  sensors  are  presented  in  Table  I.  Accord¬ 


ing  to  Table  I,  the  reconstructed  values  of  temperature  are  in 
a  good  agreement  with  the  direct  measurements:  The  dis¬ 
crepancy  is  0.1  °C  at  the  location  of  the  first  sensor  and 
0.01  °C  at  the  location  of  the  second  sensor.  These  discrep¬ 
ancies  are  smaller  than  the  estimated  errors  of  the  reconstruc¬ 
tion,  which  are  of  order  0.4  °C. 

For  the  considered  tomography  experiment,  the  tempera¬ 
ture  and  wind  velocity  fields  were  also  reconstructed  in  Ref. 
9  using  a  Simultaneous  Iterative  Reconstruction  Technique 
(SIRT)  algorithm.  SIRT  is  an  algebraic  algorithm  which  does 
not  use  a  priori  information  about  covariance  functions  of 
turbulent  fields.  It  is  worthwhile  to  compare  the  results  pre¬ 
sented  in  Ref.  9  with  those  obtained  by  TDSI.  Note  that  the 
SIRT  algorithm  used  in  Ref.  9  is  also  a  two  step  procedure. 
In  the  first  step,  the  travel  times  tf  are  split  into  tfT  and  tfv, 
where  tfT  are  affected  only  by  the  temperature  field  T, 
whereas  tfv  are  affected  only  by  the  velocity  field  V  (for 
details  see  Ref.  8).  In  the  second  step,  the  SIRT  algorithm  is 
used  to  reconstruct  the  temperature  and  velocity  fields  from 
tfT  and  tfv,  respectively.  In  Fig.  2,  the  dashed  lines  correspond 
to  the  mean  values  of  the  temperature  and  two  components 
of  wind  velocity  reconstructed  by  SIRT.  Comparison  of  these 
lines  with  the  solid  ones  (obtained  by  TDSI)  shows  that  there 
is  an  acceptable  agreement  in  reconstruction  of  70(r),  n0(r), 
and  v0(t)  by  two  algorithms.  The  root  mean  squared  discrep¬ 
ancy  between  values  obtained  by  SIRT  and  TDSI  is  0.47  °  C 
for  the  Tq,  0.16  m/s  for  the  u0,  and  0.47  m/s  for  the  u0- 
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TABLE  I.  Comparison  of  the  reconstructed  and  in  situ  measured  values  of 
temperature. 


Landscape  type 

Coordinated  (m) 

Humitter  T  (°C) 

TDSI  T  (°C) 

Bare  soil 

*=28,  >’=138 

16.24 

16.14 

Grassland 

x=  182,  >  =  143 

15.78 

15.77 

Figure  3(b)  shows  the  temperature  fluctuations  T  aver¬ 
aged  over  10  min  and  reconstructed  with  the  use  of  TDSI 
and  S1RT.  The  SIRT  reconstruction  of  T  was  taken  from  Ref. 
9  (see  Fig.  4  in  that  reference)  and  then  the  temperature  field 
was  interpolated  (2D  linear  interpolation)  to  have  the  same 
spatial  grid  as  that  used  in  TDSI.  It  follows  from  Fig.  3(b) 
that  the  reconstructed  fields  are  similar  on  large  scales  of 
turbulent  eddies.  In  the  lower  left-hand  corner  of  Fig.  3(b), 
one  can  see  a  cold  temperature  eddy  which  is  stretched  at 
about  45°  with  respect  to  the  x  axis.  Further,  there  are  two 
warm  eddies  in  the  right  upper  and  lower  corners.  However, 
there  are  differences  in  the  reconstructed  fields  in  details. 
According  to  the  TDSI  reconstruction,  the  cold  eddy  in  the 
left  lower  corner  is  more  distinctive,  narrower  and  longer 
than  that  reconstructed  by  SIRT.  The  shape  of  the  warm  eddy 
in  the  upper  right-hand  corner  is  also  more  stretched  along 
the  diagonal  and  centered  at  a  different  position.  It  is  worth¬ 
while  to  note  that,  in  TDSI,  the  deformation  of  the  turbulent 
eddies  along  the  diagonal  coincides  with  the  direction  of 
reconstructed  mean  wind  during  the  interval  of  averaging 
[black  lines  in  Fig.  3(b),  left  plot]. 

The  estimated  RMSE  for  the  TDSI  and  SIRT  reconstruc¬ 
tions  are  shown  in  Fig.  5.  Note  that  the  errors  of  the  recon¬ 


struction  by  TDSI  are  of  order  0.1  °  C,  and  are  approximately 
one-half  the  SIRT  errors. 

VI.  CONCLUSIONS 

In  this  paper,  the  TDSI  methodology  for  acoustic  tomog¬ 
raphy  of  the  atmosphere  was  further  developed.  First,  the 
formulas  for  the  spatial-temporal  covariance  functions  of 
temperature  and  wind  velocity  fluctuations  for  the  case  of 
locally  frozen  turbulence  were  derived.  These  formulas  are 
more  general  than  those  for  the  case  of  frozen  turbulence  and 
account  for  the  variance  of  velocity  fluctuations.  Second,  the 
formulas  for  the  variances  of  the  errors  in  the  input  data  in 
TDSI  were  derived. 

The  developed  TDSI  was,  for  the  first  time,  used  in  the 
reconstruction  of  temperature  and  wind  velocity  fields  in  an 
acoustic  tomography  of  the  atmosphere  experiment 
(STINHO  experiment9).  The  reconstructed  temperature  and 
wind  velocity  fields  and  the  errors  in  the  reconstruction  were 
presented.  It  was  shown  that  the  use  of  TDSI  allowed  us  to 
reliably  reconstruct  temperature  and  velocity  eddies  within 
the  tomographic  area.  The  temperature  field  reconstructed  by 
TDSI  was  compared  with  that  reconstructed  by  the  SIRT 
method.  The  fields  reconstructed  by  these  two  approaches 
are  similar  on  large  scales  while  different  in  details.  TDSI 
yields  a  more  detailed  reconstruction  with  the  errors  in  about 
one-half  those  in  SIRT. 

Note  that  the  STINHO  experiment  was  not  designed  to 
take  the  full  advantages  of  TDSI.  In  future  experiments,  it  is 
worthwhile  to  fully  exploit  these  advantages.  For  example,  it 
would  be  worthwhile  to  reduce  the  time  intervals  between 
consecutive  travel  time  measurements. 


TDSI  SIRT 
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FIG.  5.  (Color  online)  The  expected  root  mean  squared  errors  in  the  reconstructions  shown  in  Fig.  3(b). 
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APPENDIX:  SPATIAL-TEMPORAL  COVARIANCE 
FUNCTION  FOR  LOCALLY  FROZEN  TURBULENCE 


In  this  appendix,  we  derive  a  formula  for  the  spatial- 
temporal  covariance  function  of  temperature  fluctuations  for 
the  case  of  locally  frozen  turbulence.  A  starting  equation  for 
this  derivation  is  Eq.  (22).  Substituting  the  value  of  BSTT  from 
Eq.  (17)  into  Eq.  (22),  we  have 


B-nip,  T)  =  °t\  exP 


(Px~ut)2  ( Py-VT )2 


/ 

lT 


(Al) 


We  will  assume  that  the  random  3D  field  V(R ,  t)  is  normally 

2 

distributed  with  the  variance  cr~  independent  of  spatial  coor¬ 
dinates  and  time:  <x~  =(V2(R,r))  =  (w2(R,r))  +  (u2(R,f)) 
+(w2(R,f))  =  3cr^.  Its  mathematical  expectation  is  V0(f) 
=  {uQ{t)  ,u0(f)  ,w0(t)).  Note  that  different  components  of 
V(R,f)  are  uncorrelated  in  the  theory  of  homogeneous  iso¬ 
tropic  turbulence  which  is  used  in  this  paper.  Indeed,  as  tur¬ 
bulence  is  stationary,  at  any  time  t  the  cross  correlations 
between  two  different  components  of  Y,  e.g.,  u  and  v,  are 


described  byBsuv  given  by  Eq.  (19).  It  follows  from  this  equa¬ 
tion  that  the  cross  correlation  at  a  fixed  point  ( px=py=p . 
=  0)  is  equal  to  zero.  If  components  of  a  multidimensional 
normally  distributed  variable  are  uncorrelated  then  they  are 
independent.  Statistical  independence  of  different  compo¬ 
nents  of  V(R,t)  allows  one  to  express  Eq.  (Al)  in  the  fol¬ 
lowing  form: 


B-nip,  t)  =  OjIjJylz, 


(A2) 


where  /t=(exp[-(p,.-nr)2//^])i7,  and  Iy  and  Iz  are  given  by 
similar  formulas.  Let  r/={px-uT)  I  lT.  Then  Ix  can  be  written 
as  /x=(exp(-?72))^.  Note  that  a  random  function  rj  is  nor¬ 
mally  distributed  as  px,  lT,  and  r  are  nonrandom  quantities 
and  u  is  a  normally  distributed  random  function.  Therefore, 
rj  has  a  Gaussian  probability  density 


pv(w)  = - r=  exp 

(T^ZTT 


( w-(7) ))2 

- 


with  its  mean  value  ( rj)  and  variance  rr  given  by 


(v) : 


P.v~  uo(t\)T 
It 


crvr 


(A3) 


Then 


exp(-  w2)pJw)dw  = 


1 


Vl+2ol 


exp 


(v)2  \ 

1+2  </' 


(A4) 


Similar  results  can  be  obtained  for  Iy  and  l:.  Substituting 
Eq.  (A4)  and  similar  formulas  for  /,,  and  /_  into  Eq.  (A2),  we 
obtain  an  analytical  formula  for  B^p.r).  This  formula  is 
given  by  Eq.  (23). 
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